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In this paper, we will give a general introduction to the project of Ali CMB Polarization Tele- 
scope (AliCPT), which is a Sino-US joint project led by the Institute of High Energy Physics 
(IHEP) and has involved many different institutes in China. It is the first ground-based cosmic 
microwave background (CMB) polarization experiment in China and an integral part of China’s 
Gravitational Waves Program. The main scientific goal of AliCPT project is to probe the primor- 
dial gravitational waves (PGWs) originated from the very early Universe. 

The AliCPT project includes two stages. The first stage referred to as AliCPT-1, is to build a 
telescope in the Ali region of Tibet with an altitude of 5,250 meters. Once completed, it will be the 
worldwide highest ground-based CMB observatory and open a new window for probing PGWs in 
northern hemisphere. AliCPT-1 telescope is designed to have about 7,000 TES detectors at 90GHz 
and 150GHz. The second stage is to have a more sensitive telescope (AliCPT-2) with the number of 
detectors more than 20,000. 

Our simulations show that AliCPT will improve the current constraint on the tensor-to-scalar 
ratio r by one order of magnitude with 3 years’ observation. Besides the PGWs, the AliCPT will 
also enable a precise measurement on the CMB rotation angle and provide a precise test on the CPT 
symmetry. We show 3 years’ observation will improve the current limit by two order of magnitude. 


I. INTRODUCTION 


Searching for gravitational waves (GWs) has long been the cornerstone of cosmology and astrophysics since Einstein 
proposed the General Relativity (GR) in early 20th century. GWs are thought to be the last piece of the theoretical 
predictions of GR. With long-lasting efforts, LIGO (Laser Interferometer Gravitational-Wave Observatory) collaboration 
in 2016 announced the first detection of GWs with the signals coming from two merging black holes with mass of tens 
of solar masses [I]. Since then, LIGO and Virgo have announced other three events of black hole GWs [2H4]. These 
achievements make the GWs study to enter a new era and were awarded the Nobel Prize in Physics this year. 

Different from the GWs detected by LIGO and Virgo, the PGWs arise from quantum fluctuations and carry 
important information about the very early Universe, for example, the physics of inflation, bouncing and emergent 
Universe. So far the most effective way to probe PGWs is to measure the B-mode polarization of CMB. 

The CMB photons are relics left after the Big Bang. Its first detection half a century ago pioneered the study of 
cosmology. For the recent 20 to 30 years, CMB observations have developed rapidly, leading us into the precision 
cosmology era. However, the CMB B-mode polarization induced by the tensor fluctuations generated in the early 
Universe, i.e., the PGWs, still have not been detected conclusively!. This has become a key scientific goal of CMB 
observations in recent years. In addition, CMB B-modes provide us an important test on fundamental physics, such as 
the CPT symmetry. 


* Invited Paper for National Science Review 
1 BICEP2 collaboration[5] announced in 2014 the detection of B-mode, however it turns out to be the dust emissions , not the signal of 


PGWs[6]. 


At present, major ground-based CMB experiments are in southern hemisphere, for example, the Atacama Cosmology 
Telescope (ACT) and POLARBEAR/Simons Array in Chile, and the South Pole Telescope (SPT) and BICEP at the 
South Pole. High precision experiments in the northern hemisphere are critically needed to achieve a full sky coverage. 

In 2014, IHEP cosmology team proposed a CMB experiment in Ali of Tibet, aiming to search for the PGWs in 
northern hemisphere. In this paper we provide a general introduction to the AliCPT project. In section|II| we introduce 
the atmospheric conditions, sky coverage and infrastructure of AliCPT. In section [II] we present the main scientific 
goals of AliCPT. Section [IV] is our summary. 


II. OVERVIEW OF THE ALICPT SITE 


Ali prefecture is in the west of the Tibetan Plateau with vast highlands and the site of AliCPT is located on top of 
a peak in the northwest of Gangdise mountain range with many peaks higher than 5,000m and even 6,000m. The 
Xinjiang-Tibet Highway (officially called the China National Highway 219) as well as the Ngari Gunsa Airport are not 
far from the site. 

In this section, we describe the atmospheric conditions, sky coverage and the infrastructure of AliCPT site, which 
are crucial to ground-based CMB experiments. 


e Atmospheric conditions: For ground-based CMB telescopes, the atmosphere is inescapable. The absorption 
and emission at millimeter/sub-millimeter band by the air molecules reduce the significance of signals. Among 
all components of the air, water vapor plays a crucial role due to its strong absorption/emission and heavy time 
variation. Usually, the Precipitable Water Vapor (PWV) is used as a conventional parameter to characterize the 
amount of water vapor, which is defined as the overall depth of water in a column of the atmosphere above the 
ground. CMB signals, especially the polarization signals, are extremely weak, so observation requires the air of 
the site to be thin, dry and stable. In the upper panel of Figure [I] we show the global distribution of the mean 
values of PWV over the past 6 years (2011.7-2017.7). It’s obvious that only four regions have the lowest PWV 
on earth, including the Antarctic, Atacama Desert, Greenland and the high Tibetan Plateau. 


AliCPT site is located at a 5,250m high peak of the Gangdise Mountain. The Himalayas is to its southwest 
and runs from northwest to southeast, separating Ali area from the Indian subcontinent as well as the Indian 
Ocean, which is shown in the lower panel of F igure [I] So the wet air from the Indian Ocean is largely reduced. 
This makes air of AliCPT site to be thin and dry enough around Winter. In [7], we quantitatively analyzed the 
atmospheric conditions of the site using radiosondes data from the local weather station as well as MERRA-2 
reanalysis data from NASA/GMAO. The results show that PWV of Ali have a very strong seasonal variation and 
the median PWV of the observing season (October to March) is about 1mm (1.07mm for MERRA-2, 0.92mm 
for radiosondes), which is excellent to observation at 90/150GHz. We show cumulative fraction of time versus 
PWV in Figure [2| References [8] [9] have also done the atmospherical conditions evaluation for Ali region. 


e Sky coverage: The AliCPT site is located at geographical coordination (80°01’E , 32°19’N). With the 
rotation of the earth and the mid-latitude location, AliCPT is able to cover the whole northern sky as well as 
the low latitude part of the southern sky, and the overall observable fraction is about 70 percent. We show this 
in Figure [3] as the region above the black dashed line. In our calculation we choose the instrumental parameters 
to be a 45° lowest elevation for the mount and a 30° field of view (FOV). The overlap in the low latitude region 
of observable sky of Ali and Atacama makes it convenient to do cross-check and cross-correlation studies. Known 
as the northern hole, the lowest foreground contaminated region in the northern galactic hemisphere is also 
within the observable sky of AliCPT, which is extremely important for a CMB B-mode polarization and PGWs 
aimed project. We show the target fields of AliCPT in Figure [3] TN1 and TN2 within the black solid lines are 
target fields in the northern galactic hemisphere and TS in the southern galactic hemisphere. TND with the 
lowest dust intensity is chosen for deeper survey. 


So the sky coverage (both observable and low foreground contamination sky) of AliCPT is complement to that of 
experiments in Antarctic and Atacama. Together with the southern cleanest sky region that covered by southern 
projects such as BICEP and Simons Array, AliCPT will increase the chance to find the B mode and PGWs. For 
more details about sky coverage, we refer to [7] . 


Infrastructure: On the transportation, the National Highway 219 is right next to the AliCPT site. Moreover, 
the Ngari Gunsa Airport is only about half an hour driving distance to the site and has a daily commercial flight 
to Lhasa, the capital of Tibet. The largest settlement of Ali area, Shiquanhe town is also very close to the site 
and it takes only 30 minutes driving to get there. 
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FIG. 1: Global distribution of mean PWV over 6 years (2011.7-2017.7), obtained with MERRA-2 data (upper). The location of 
AliCPT site and the terrain around. The wet air from the Indian Ocean is hugely reduced by the Himalaya Mountain(lower). 


From the aerial view around Ali astronomical observatory at Al point (left panel of Figure|4), we can see that, 
AliCPT site at B1 is not far from Al point, the distance between them is only about 1km. The concrete road 
from A1 to B1 is already under construction and will be finished soon. In Ali astronomical observatory at A1 
point, city grid electric power and the network infrastructure for data transmission are ready now. Some optical 
telescopes for astrophysics have been set up. Started in March 2017, the site construction of AliCPT is already 
ongoing and will be finished by the end of this year. After the commissioning in 2019, the observation is expected 
to start in 2020. 


In short, AliCPT opens a new window in the northern hemisphere to detect CMB B-mode polarization and probe 
the PGWs. 
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FIG. 2: The cumulative distribution for sites at 5,250m and 6,000m over observing season (left) and comparison between 
MERRA-2 and radiosondes datasets (right). Results are taken from [7]. 
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FIG. 3: Sky coverage of AliCPT. The observable sky is above the black dashed line, the target field is shown as TN1 and TN2 
in the northern galactic hemisphere and TS in the southern galactic hemisphere. TND is for deeper survey. 


III. ALICPT AND ITS SCIENTIFIC GOALS 


In this section, we firstly introduce the main contents of AliCPT project, then perform the simulations on scientific 
goals. AliCPT project consists of two stages. The first stage is to develop and deploy a CMB polarization telescope 
at 5,250 meters, called AliCPT-1. AliCPT-1 telescope is a dichroic refractor of aperture 70cm covering 90/150GHz, 
with a three-axis driving mount scanning at speed of 5°/s in azimuth. AliCPT-1 adopts Transition Edge Sensor(TES) 
bolometers used widely in current CMB polarization experiments{10], and takes Superconducting Quantum Interference 
Devices (SQUIDs) as cryogenic readout. Sensors and their readout will be packaged into highly integrated modules, 
and each module contains 1,704 TES sensors. Four modules will be installed in AliCPT-1 telescope by the end of 2019 
and the number of detectors will reach 6,816. 
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FIG. 4: From left to right: aerial view of the region around Ali Astronomical Observatory, picture for laying the foundation for 
AliCPT at B1 point, the current infrastructure for AliCPT. 


TABLE I: Instrumental parameters for AliCPT-1 and AliCPT-2. Schedule for the number of modules and detectors installed. 
NET is short for Noise Equivalent Temperature, fsky represents sky coverage. 


Year 2019 (aticpr-1) | 2020(alicpT-1+AlicPT-2) | 2021 (A1iCPT-1 + AlicpT-2)|2022(alicpT-1 + AliCPT-2) 
NET(uK vs) 350 350 350 350 

Nmoa 4 4+4 4+8 4+12 

Naet 6,816 13,632 20,448 27,264 

fy 10% 10% 10% 10% 
Bands(GHz) | 90 & 150 90 & 150 90 & 150 90 & 150 


In the second stage, we will have a more sensitive telescope (AliCPT-2) with 12 modules and detectors more than 
20,000. The construction of AliCPT-2 will start in 2020. In each year, we will install four modules and finish the 12 
modules by the end of 2022. In Table I, we summarize the basic instrumental parameters used for the simulations and 
the corresponding schedule of AliCPT-1 and AliCPT-2. 

The main scientific goals of AliCPT include: 


e Start the large coverage surveys in northern hemisphere and search for regions with low foreground contamination; 
e Target at the cleanest sky region and detect PGWs in the northern hemisphere; 

e Measure the CMB rotation angle in high precision and test the CPT symmetry; 

e Study the hemispherical asymmetry in combination with experiments in southern hemisphere; 

e Measure E-mode polarization in high precision and study the effects in cosmology; 

e Study cross-correlation of CMB polarization with the large scale structures (LSS). 


In the following, we provide the simulation on the sensitivity of the measurement on r and the CMB rotation angle. 


A. Sensitivity on r and its implication on the Early Universe Physics 


The leading scenario of the early Universe is the inflationary cosmology, which describes an accelerating expanding 
phase occurred before radiation epoch. Inflation resolves several conceptual issues of the Big Bang theory including 
the flatness, monopole, and horizon problems{Ii]. Moreover, inflation explains the origin of primordial perturbations, 
with a mechanism that the quantum fluctuations of inflaton field were stretched to be classical perturbations by 
exponential expansion of space. The primordial perturbations have three types: scalar, vector and tensor. Scalar 
modes will eventually seed the CMB temperature anisotropies and also lead to the formation of large-scale structures 
in the Universe. Tensor modes, which are dubbed as PGWs, will introduce CMB B-mode polarization, which is the 
target signal in AliCPT observations. Conventionally, we often use r, ng and some other parameters to describe tensor 
perturbations, where r means the ratio of the amplitudes of power spectra of primordial tensor Ar and scalar modes 
Ag, n represents the logarithmic derivative via 


Ar _ dlnP, 
~ dink ’ 


(1) 
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FIG. 5: Comparison of theoretical model predictions with the CMB experiments through the n; - r diagram. The green line and 
red line give the 20 sensitivity reached by AliCPT until 2022 and 2025. 


where P, is the power spectrum of tensor mode, k is the wave number. 

Since inflationary cosmology still suffer from the initial cosmic singularity issue[12] [I3], alternative theoretical 
approaches have been proposed, including bounce cosmology [I4], cyclic Universe [15], and emergent Universe [I6]. 
Similar to inflation, these theories can also generate primordial tensor perturbations. However, PGWs generated 
from different theories have different characteristics, which will be reflected on the shape of CMB BB spectrum. The 
comparison of typical model predictions for the r — n, relationship with the present and forthcoming CMB experiments 
is shown in Figure [5|17]. We can see that high precision measurements of PGWs are crucial in testing models of early 
Universe. 

We have performed simulations to forecast the constraining capability of the AliCPT survey. To obtain the constraint 
on r derived from the AliCPT observations, we adopt the Fisher matrix approach [18], which is an efficient way to 
forecast the constraints on cosmological parameters given the specification of the instruments. 

Generally, the likelihood function of a multivariate Gaussian-distributed data vector d can be expressed as, 


1 


l 4 
L= CO) exp ( 5d [C(A)} a) i (2) 


where @ is the parameter vector, and C is the covariance matrix which is a function of 0 in general. The Fisher matrix, 
Fij, which is the second partial derivative of the likelihood function with respect to parameters 6; and 0; evaluated 
at the fiducial model, approximates the Hessian matrix, and the inverse of the diagonal terms (F~!),;; provides an 


estimate of the lower limit of the variance for parameter 0;, ie., A@; Z (F~!)}/? [I8]. 
For the case of CMB, d = {X7,,,...}, where X € {T, E, B}, snd v runs over all available frequency bands. Each 
X7,,, consists of three components, namely, the lensed CMB, the foreground emission and the instrumental noise. The 


Fisher matrix for the CMB observables is, 


20 1 Cpt OO 
Fig = >> fay Tr [or 99, Ce E 30, (3) 
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where £ denotes the order of the multipole, fsky is the sky coverage, and Cy is the covariance matrix in harmonic space 
which can be expressed as, 
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FIG. 6: Left panel: Sensitivity of r without delensing. Right panel: Sensitivity of r with completely delensing. The horizontal 
black dashed line is the current limit r < 0.07. 


where Ce, Fe, Neg correspond to spectra of the lensed CMB, foreground and instrumental noise respectively. We 
have assumed that observed temperature and polarization are statistically isotropic, so that terms of different @’s are 
independent. 

We use CAMB [I9] to calculate the CMB spectra Ce, and choose a fiducial cosmology which is consistent with the 
on Planck 2015 results [20]. Theoretically, the B-mode signal from lensing will exceed that from the primordial signal 
aa when on recombination bump (£ ~ 100) r falls below 0.01. Therefore, to achieve a high precision detection of the 
primordial B-mode signal, we shall perform a delensing procedure to remove the lensing effect from data. We consider 
two cases in our simulations for a comparison, namely, the cases of totally-delensed and undelensed. During the actual 
operation, we will use the data of AliCPT and other small scale CMB surveys at high £ to reconstruct lensing spectra. 
Other approaches, e.g., cross-correlating CMB with large scale surveys, will also be considered. 

Based on results of recent experiments, foreground emission Fẹ is thought to dominate over all frequency bands 
and all scales. This kind of contamination can be removed through multi-frequency surveys, since the frequency 
spectrum distributions of CMB and other components are different. Component separation is then defined as a 
method estimating the amount of each emission component. Even such a separation has been applied, some extent 
of residual foreground still remains in data. The AliCPT telescope plans to scan the cleanest sky regions at 90/150 
GHz, in order to suppress residual foreground as much as possible. In our Fisher forecasting, we consider two kinds of 
components, namely, the synchrotron and thermal dust, which contribute to the dominate part of polarized foreground 
emission. Besides r, we take six additional parameters into free parameter space, including dust, Async corresponding 
to amplitudes of dust and synchrotron at @ = 80 with pivot frequencies of 353 GHz and 23 GHz respectively, @aust and 
Qsync Corresponding to spectral index in harmonic space, faust and sync corresponding to spectral index in frequency 
space. And no conversion between these two components, dust and synchrotron is assumed in our study. 

In our simulations, we assume that the AliCPT instrumental noise of temperature and polarization are uncorrelated, 
and an isotropic Gaussian random noise for simplicity. The temperature spectrum is given by Ne = w~'B? where By 
is the harmonic transform of the Gaussian beam, reads, By = exp(—£(¢ + 1)62wyy/8 log 2). The weight is then, 


-1 _ 4T fsky NET? 
tobsNdet 


w (5) 
where NET means the Noise Equivalent Temperature which is closely related with the detector performance, instrument 
design and atmospheric condition, tops is the effective observation time (from October to March, 10h per day) and Naet 
is the number of detectors. Note for polarization, NET should be multiplied by a factor of v2 since each polarized 
signal needs two orthogonal linear polarized detectors. In Table |I| we show the instrumental parameters for our 
calculation. 

We plot our results in Figure [6] Simulated sensitivity of r under undelensed and totally-delensed assumptions are 
showed in the left and right panels. Cases related to 1, 5, and 10 percent of residual foreground are considered in both 
panels. Dashed black lines represent current representative constraint on r, r < 0.07 at 20, obtained from the joint 


TABLE II: Measurements on CMB rotation angle since 2006. 


Data a +o 4 Sus 
1 |}WMAP3+BOOMERANG —6° + 4° 
2 WMAP3 BI] —2.5° + 3.0° 
3 WMAP5 -1.7° +2.1° 
4 WMAP7 —1.1° + 1.4° +1.5° 
5 WMAP9 [Bd] —0.36° + 1.24° + 1.5° 
6 QUaD —0.56° + 0.82° + 0.5° 
7 BICEP1 —2.6° + 1.02° 
8 BICEP1 —2.77° + 0.86° + 1.3° 
9 POLARBEAR [38] —1.08° + 0.20° + 0.5° 
10 ACTPol —0.2° +0.5° 
11 Planck 2015 0.35° + 0.05° + 0.28° 


analyses of BICEP2/Keck Array and Planck[2I]. As one can see from the plots, even in the undelensed case with 
10 percent residual foreground, after 3 years’ survey of AliCPT, the sensitivity on r will reach o, = 0.007, which is 
nearly one order of magnitude stronger than current constraint. By the end of 2025, after 3 years’ observation with 16 
detector modules (around 27000 TES detectors), the constraint on r in delensed case will reach a level of 0.003. Not 
surprisingly, less residual foreground leads to a better result on r. If we can suppress foreground down to the one 
percent level, and apply a delensing procedure, we will obtain a, ~ 0.001. Much stronger limit of r can help us with 
testing a large number of cosmological models of the early Universe. Due to its high precision, the AliCPT project will 
provide us a new vision of the early Universe. 


B. Sensitivity on the CMB Polarization Rotation Angle and its implication on the CPT Test 


Testing CPT symmetry, the combination of charge conjugation (C), parity reflection (P) and time reversal (T) is 
important to cosmology and particle physics. Any violation, if found, would be a powerful and important clue for 
new physics beyond the standard model. So far, the laboratory experiments are consistent with a null result for the 
CPT violation. However, these tests may not be applied to physical processes in the early Universe at extremely high 
energy scales. In fact, there are motivations to speculate on the CPT violation in cosmology. Firstly, the expanding 
Universe has a preferred temporal direction, which provides a natural frame to break the Lorentz and CPT symmetries. 
Secondly, the observed baryon and anti-baryon asymmetry in the Universe may indicate a dynamical CPT violation 


[22 23]. 


To study the cosmological CPT violation in CMB, we consider the following effective Lagrangian, 
Los = pA F" , (6) 


where the external field p, is a constant vector[24], or p, ~ ô „$ with @ being the dark energy scalar in quintessential 
baryo/leptogensis[23] [25], or p, ~ „R with R the Ricci Scalar in gravitational baryo/leptogensis|[26] [27], Fw — 
(1/2)e"”"? Foo is the dual tensor of the electromagnetic tensor Fy. With the Chern-Simons term in (6), the polarization 
directions of photons get rotations for CMB, this will convert part of E-mode polarization to the B-mode and change 
the power spectra of polarization fields [28H30]. In Table [II] we summarized the constraints on the rotation angle from 
various experiments. The current limit on the rotation angle a is about 1°. 

We have performed simulations to forecast the sensitivity on the measurements of the rotation angle with AliCPTs. 
For the forecast, we employ the so-called D-estimators[41]| defined as follows, 


DIBebs — CTB.OPS 695/98) — CTE sin(28); 
il 
DEB. obs = CBBiobs daij = (0p _ CBB.obs) sin(48), (7) 


where £ is the unbiased estimator for the cosmic rotation angle a. 
We modified the generic Monte Carlo Markov Chain (MCMC) sampler provided in the CosmoMC package|42] and 
performed the calculation with instrumental properties of AliCPT. We present our results in Figure [7] where the blue 
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FIG. 7: Forecast of average polarization rotation angle. The blue dashed line is obtained from TB estimator and the black solid 
line is from TB+EB estimator. 


dashed and black solid lines represent the constraints with TB and EB + TB estimator respectively. Our results show 


that, after 3 years’ observation, the AliCPTs are capable of providing a stringent constraint on the average rotation 
angle of o(a) ~ 0.01°. 


IV. SUMMARY 


In this paper, we have provided an introduction to the AliCPT project and its scientific goals. After its completion, 
together with Chile and the South Pole, AliCPT observatory will be one of the major CMB polarization observation 
stations in the world, realizing a full sky coverage in the search of the PGWs. 
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